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I.  SUMMARY 

Diamond  is  a  material,  exhibiting  unique  mechanical,  thermal  and  electrical  properties,  as  well  as  high 
electron  and  hole  mobility  [1],  promoting  its  high  performance  in  microelectronic  devices.  In  this  study 
we  have  investigated  the  microscopic  processes  involved  in  laser  irradiation  of  diamond  and  then  we  have 
extended  it  for  three  different  potential  applications. 

Photoionization  of  diamond  below  and  close  to  the  graphitization  threshold  as  in  the  experimental  pro¬ 
cedure  used  to  write  graphitic  electrodes  inside  diamond  bulk  by  laser  irradiation  [2—4]  for  the  fabrication 
of  three-dimensional  (3D)  devices  was  theoretically  described  and  compared  to  laser-induced  current  in 
monocrystalline  diamond  detector  measured  by  Transient  Currents  Technique  (TCT)  in  a  wide  range  of 
laser  intensities  and  at  different  bias  voltages. 

In  connection  to  the  fact  that  diamond  is  a  basic  ingredient  in  modern  nanophotonics  [5,  6]  and  due  to 
its  high  refractive  index  in  UV-VIS  range,  it  is  prospective  material  for  all-dielectric  [7]  and  hybrid  metal- 
dielectric  nano-photonic  devices  and  circuits  [8-10]  we  next  have  studied  the  process  in  which  despite  its 
dielectric  character,  similarly  to  silicon  it  can  be  promptly  turned  by  intense  ultrashort  laser  pulses  into 
short-lived  plasmonic  state,  becoming  so-called  ’’virtual  plasmonic  material”,  supporting  photoexcitation 
and  propagation  of  surface  plasmon-polaritons  (SPPs).  The  potential  applications  are  in  ultrafast  optical 
switching,  spatial  phase  modulation  and  saturable  absorption  [1 1-15]. 

The  feasibility  for  femto  and  attosecond  pulse  generation  from  diamond  subjected  to  intense  fs  laser 
pulse  with  mid-infrared  wavelength  was  explored  by  further  investigating  the  non-perturbative  scaling  of 
the  generated  harmonics  with  the  laser  field  and  their  bandstructure  dependence  since  recent  experimental 
and  theoretical  investigations  in  the  literature  show  that  high  harmonic  generation  HHG  in  wide  band-gap 
semiconductors  (insulators)  may  lead  to  brighter  sources  due  to  the  high  laser  induced  carrier  density  and 
ability  to  withstand  the  driving  laser  fields  necessary  for  inducing  strong  nonlinear  and  non  perturbative 
effects  leading  to  HHG  [16-  0] 


II.  PHOTOIONIZATION  OF  MONCRYSTALINE  CVD  DIAMOND  IRRADIATED  WITH 
ULTRASHORT  INTENSE  LASER  PULSE 

A.  Introduction 

Writing  graphitic  electrodes  inside  the  diamond  bulk  by  laser  irradiation  is  achieved  in  the  femtosecond 
range,  from  400  fs  to  30  fs  with  no  damage  of  the  diamond  lattice  that  has  been  reported  to  occur  in  the 
nanosecond  range[21].  This  technique  is  used  for  producing  radiation  tolerant  three-dimensional  diamond 
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detectors  in  high  energy  physics  and  small  field  dosimeters  for  radiotherapy  in  which  a  shorter  interelectrode 
distance  results  in  a  lower  voltage  bias  required  to  collect  the  charge  produced  in  the  diamond  thickness. 
Because  of  the  reduced  path  of  the  charge  earners,  the  charge  loss  is  minimized  and  the  Charge  Collection 
Efficiency  (CCE)  increased,  as  a  consequence,  these  detectors  can  be  can  be  the  radiation  hardest  sensors 
ever  implemented  [22],  Only  one  theoretical  model  on  diamond  graphitization  has  been  proposed  so  far, 
involving  single  photon  excitation  above  bandgap[23,  24],  which  is  an  experimental  condition  distant  from 
the  cited  experimental  conditions  and  applications,  because  irradiation  below  bandgap  is  needed  to  produce 
electrodes  inside  the  diamond  bulk. 

A  deep  understanding  of  the  physics  of  the  processes  is  required  to  predictively  tune  the  laser  parameters 
(energy  per  pulse,  pulse  width,  repetition  rate),  in  order  to  obtain  the  sought  properties  of  the  electrodes 
inside  the  diamond  bulk  such  as  resistivity  and  charge  collection  efficiency.  These  are  the  processes  of 
excitation  and  recombination  of  carriers,  release  of  energy  to  the  lattice  and  modification  of  the  chemical 
bond.  It  is  also  critical  to  take  into  account  the  different  crystalline  qualities  of  the  diamond  samples 
available,  which  can  be  determined  by  then-  original  optical  and  electronic  properties. 

We  used  a  quantum  theoretical  model  to  simulate  multi-photon  excitation  across  the  diamond  bandgap. 
At  the  highest  energies  the  change  in  optical  parameters  of  the  diamond  caused  by  the  high  concentration  of 
earners  was  taken  into  account.  Non-linear  absoiption  causes  a  relevant  amount  of  energy  to  be  deposited 
prior  to  the  focus  as  observed  in  other  materials  [25]  .  As  a  consequence  the  energy  density  required  to  the 
phase  transition  in  the  pure  diamond  bulk  cannot  be  reached.  We  have  shown  that  the  theoretical  model 
supplemented  with  high  energy  correction  is  in  agreement  with  the  experimental  data  over  a  very  wide 
energy  and  corresponding  carrier  density  range,  resulting  in  the  first  quantitative  validation  of  a  theoret¬ 
ical  model  of  multi-photon  excitation  of  diamond  to  our  knowledge.  This  result  opens  the  way  to  a  full 
understanding  of  laser  graphitization  by  femtosecond  pulses,  a  subject  relevant  also  beyond  the  suggested 
applications.  Diamond  is  a  cardinal  material  with  ever  growing  experimental  applications  and  the  process 
of  laser  graphitization  is  generally  important,  even  when  it  is  an  unwanted  side  effect. 


B.  Methods,  Assumptions  and  Procedures 

For  the  relatively  short  duration  of  the  laser  pulse  (30fs)  considered  in  this  work  we  assume  the  nuclei 
remain  fixed  at  the  lattice  points  during  the  irradiation,  such  that  energy  is  absorbed  from  the  laser  field  via 
creation  of  electron-hole  pairs.  We  also  neglect  the  contribution  of  impact  ionization  to  the  population  of 
the  conduction  band.  It  has  been  predicted  that  avalanche  generation  in  dielectrics  is  less  important  than 
strong-field  ionization  for  pulse  durations  below  100  fs[26,  27],  However,  there  is  experimental  evidence 
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that  optical  breakdown  in  fused  silica  is  dominated  by  avalanche  ionization  down  to  the  10  fs  regime[28].  In 
particular  the  observed  multiphoton  ionization  rates  were  found  substantially  lower  than  those  predicted  by 
Keldysh  theory,  thus  identifying  non-perturbative  mechanism  of  earner  generation  in  bulk  silica.  Moreover, 
in  a  recent  work  on  non  linear  absorption  of  femtosecond  laser  pulses  in  silica  [29]  transmission  curves  at 
different  pulse  durations  have  been  fitted  to  a  model  taking  into  account  multiphoton  absorption  and  impact 
ionization  through  a  modified  rate  equation.  The  results  were  consistent  with  a  photon-assisted  avalanche 
ionization  rate  that  increases  with  decreasing  pulse  length  down  to  a  10  fs  time  scale.  Nevertheless  the 
physical  conditions  under  which  the  prominence  of  impact  ionization  occurs  in  silica  will  be  shown  to  be 
very  different  from  those  encountered  in  diamond  at  the  same  laser  intensities.  In  fact,  even  if  the  kinetic 
energy  distribution  of  the  photoexcited  carriers  extends  above  the  required  threshold,  the  process  of  impact 
ionization  is  strongly  suppressed  by  the  high  population  of  the  conduction  band  low  energy  states  due  to 
Pauli  blocking  and  neglecting  the  impact  ionization  in  our  model  still  provides  good  quantitative  agreement 
with  the  measured  photoionization  yield  in  diamond. 

More  recently,  it  was  experimentally  found  that  high-harmonic  generation  in  bulk  solids  [30]  is  char¬ 
acteristic  of  a  non-perturbative  quantum  interference  process  that  involves  electrons  from  multiple  valence 
bands.  Thus  in  order  to  clearly  distinguish  field-assisted  collisional  effects  from  coherent  process  in  the 
laser-irradiated  diamond,  we  investigate  the  quantum  dynamics  of  the  microscopic  polarization  and  charge 
carrier  distribution  in  one -particle  framework  by  incorporating  realistic  substrate  band-structure. 

In  long  wavelength  approximation  we  represent  the  light  pulse  by  a  spatially  uniform  time-dependent 
electric  field.  The  coupling  of  electrons  to  the  laser  field  is  modeled  in  single-particle  approximation  based 
on  the  Hamiltonian 

H(t)  =  ^(p  +  A(t))2  +  V(r)  (1) 

including  the  periodic  ion-lattice  potential 

V (r)  =  V (G)  cos(G  •  r )eiG'r,  (2) 

G 

where  2 r  =  ao(l/4, 1/4, 1/4)  is  a  relative  vector  connecting  two  carbon  atoms  in  a  crystal  unit  cell, 
ao  =  3.57.4  is  the  bulk  lattice  constant  for  diamond  and  G  labels  the  reciprocal  lattice  wave-vectors.  We 
apply  the  empirical  pseudopotential  method  to  describe  the  atomic  form  factors  V(G)  [31,  32].  Velocity 
gauge  implied  by  Eq.(l)  is  used  throughout  our  calculation  [33]  and  the  time-dependent  electronic  wave- 
function  is  expanded  over  Bloch  states  of  definite  crystal  momentum  k 

V’O; f)  =  ank(f)eik'runk( r),  (3) 

nk 
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where  un k  is  the  lattice-periodic  part  of  the  Bloch  wave-function 

w„k(r)  =  ^2  e*G'rcn,G+k,  (4) 

G 

expressed  in  terms  of  Fourier  coefficients  satisfying 

-(G  +  k)2Cnic+k  +  2L  ^G-G,cn,G'+k  =  enkCn,G+k,  (5) 

G' 

with  eigen-energies  cr,k  specifying  the  static  substrate  band  structure.  The  time-evolution  of  the  Fourier 
amplitudes 

idtan kk)  (t)  =  £nk«,(!'kk)  (t)  +  P W  (k)a;™k)  (t) ,  (6) 

n' 

is  subject  to  initial  conditions  specified  long  before  the  staid  of  the  pulse  a;^k)  (t  — >  —  oo)  =  5nrn  and  the 
matrix  representation  of  the  momentum  operator  is 

Pnn'(k)  =  kdnn’  -  i- —  [  d3ru*k(r)Vr«n'k(r),  (7) 

“«.c.  Jau. c. 

with  off-diagonal  matrix  elements  describing  interband  transitions,  flu  c.  =  «jj/4  is  the  volume  of  the  unit 
cell.  The  occupation  numbers  of  single-particle  states  at  time  t  arc 

M 

fnk{t)  =  laikk)(*)|2>  (§) 

m=  1 

where  M  is  the  total  number  of  in  initially  occupied  valence  bands.  The  total  number  of  electrons  that  have 
been  photo-excited  into  the  conduction  band  is  given  by  the  Brillouin  zone  (BZ)  integral 

°SL  r 

ne{t)  =J2  (9) 

n=M+l  JbZ  ^ 

Charge  neutrality  implies  that  the  total  number  of  valence  band  holes  is  equal  the  total  number  of  conduction 
electrons,  i.e.  rih  =  ne. 


C.  Results  and  Discussion 

We  evaluate  the  substrate  band  structure  in  a  plane-wave  basis  with  kinetic -energy  cut-off  35  Ry.  The 
pseudopotential  form  factors  V(G)  used  for  diamond  (in  Rydbergs)  are  V(G2  =  3)  =  —0.625,  V(G2  = 
8)  =  0.051  and  V(G2  =  11)  =  0.206,  here  the  wave-number  G  is  given  in  units  of  2n/ao.  The  time- 
dependent  electron  density  is  represented  in  a  basis  of  static  Bloch  orbitals  including  4  valence  and  16 
conduction  bands,  and  the  time-integration  of  the  Schrodinger  equation  is  performed  by  applying  the  Crank- 
Nicolson  method  for  small  equidistant  time  steps  St  «  1  attosecond.  The  BZ  integration  is  carried  out  by 
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a  simple  Monte  Carlo  method  using  1500  randomly  generated  k-points  in  a  cube  of  edge  length  (47r/ao). 
The  time -profile  of  the  laser  field  is  modeled  by  a  Gaussian  function 

A  (f)  =  eAo  sin(wif)e_ln('4')t2/Ti,  (10) 

where  e  is  a  unit  vector  pointing  in  the  direction  of  the  electric  field,  the  photon  energy  related  to  the  driving 
frequency  is  Hujl  =  1-55  eV,  the  pulse  duration  is  tl  =  30  fs  and  the  field  strength  A$  is  given  by  the  peak 
intensity  of  the  laser  field  with  I  =  cu^AfJ (87r)  where  c  is  the  speed  of  light  in  vacuum. 

In  Fig.l,  we  show  the  static  band  structure  for  crystal  momentum  k  changing  along  the  TX  and  TL 
directions  in  the  bulk  BZ.  The  empirical  pseudopotential  model  reproduces  quantitatively  the  principal 
energy  gaps  of  this  band  structure.  The  location  of  the  conduction  band  minimum  at  k  «  (0.8, 0, 0)2-7r/ao 
corresponds  to  an  indirect  gap  of  5.42  eV,  that  is  in  very  good  quantitative  agreement  with  the  experimental 
result  (Cgap  =  5.45  eV).  Since  in  long-wavelength  approximation  the  laser  cannot  induce  a  momentum 
changing  transition,  the  minimal  energy  required  to  surpass  the  energy  gap  at  the  T  point  is  7  eV.  For 
relatively  low  intensity  of  the  laser  field  I  <C  100  TW/cm2,  excitation  of  valence  electrons  across  this  gap 
should  occur  via  5  photon  absorption  according  to  the  Keldysh  theory [34], 

In  Fig. 2  we  plot  the  time  evolution  of  the  electron  density  ne(t)  for  linearly  polarized  laser  field  directed 
parallel  to  the  bonds  between  carbon  atoms  and  for  three  different  intensities  I  =  1.5  TW/cm2,  I  =  10 
TW/cm2  and  I  =  50  TW/cm2.  The  total  number  of  electrons  excited  to  the  conduction  band  increases  on 
the  rising  part  of  the  pulse.  At  a  later  time,  the  charge  density  oscillations  follow  the  time  periodicity  of 
the  incident  radiation.  These  transient  charge  oscillations  are  exclusively  due  to  the  time-dependent  shift 
of  the  momenta  of  valence  electrons  k  — >•  k  +  A(t),  which  adjust  instantaneously  (adiabatically)  to  the 
slow  variation  of  the  vector  potential.  Due  to  this  velocity  boost,  the  overlap  among  shifted  valence  and 
undistorted  conduction  bands  (ck|nk(t))  oscillates  with  the  periodicity  of  the  laser  field.  During  one  half¬ 
cycle  of  the  pulse,  some  fraction  of  electrons  is  transferred  back  into  the  valence  band,  and  deviation  from 
the  adiabatic  time-evolution  occurs  at  the  extrema  of  the  electric  field.  The  photo-ionization  yield  stabilizes 
in  the  wake  of  the  pulse.  Photoionization  is  less  likely  for  lower  field  intensities  ( I  <  10  TW/cm2),  since 
valence  electrons  do  not  gain  sufficient  energy  from  the  field  to  be  excited  into  the  conduction  band.  That 
is  because  of  the  wide  bandgaps  of  diamond,  which  favor  adiabatic  time  evolution  resulting  in  relatively 
small  final  electron  yields  ne  ~  1018  cm-3.  Non-adiabatic  effects  of  electron-hole  pair  creation  depend 
very  sensitively  on  the  laser  intensity.  When  the  intensity  is  raised  above  I  >  50  TW/cm2  the  occupation 
of  the  conduction  band  becomes  prominent  resulting  in  high  photoionization  yields  with  ne  ~  1 02 1  cniL 3 
after  the  conclusion  of  the  pulse. 

In  Fig. 3  we  plot  the  time-dependent  macroscopic  polarization  density  P(t)  =  f1  ,x  dt'J(t'),  which  gives 
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the  charge  transferred  into  the  conduction  band  per  unit  area  at  time  t.  The  macroscopic  current  density  J(t) 
in  the  direction  of  the  electric  field  is  given  by  an  average  over  the  volume  of  the  unit  cell 

J(t)  =  -^—[  d3r  j(r,  f)  •  e,  (11) 

where  the  space-  and  time-dependent  electron  current  density 

jO,  *)  =  ^2  \  [Ck(r.  t)(-iVr  +  A(t))V>mk(r,  t)  +  c.c]  ,  (12) 

mk 

is  given  by  the  sum  over  initially  occupied  valence  band  orbitals.  For  relatively  low  intensity  of  the  laser 
pulse  /  ~  2  TW/cm2,  the  polarization  follows  adiabatically  the  applied  electric  field.  For  the  increased 
pulse  intensity  I  =  10  TW/cm2,  the  dielectric  response  of  the  crystal  remains  predominantly  linear  with 
no  noticeable  sustained  current  of  delocalized  electrons  after  the  end  of  the  pulse.  When  the  intensity  is 
increased  to  I  =  50  TW/cm2,  the  polarization  density  exhibits  almost  linear  time  dependence  in  the  wake 
of  the  pulse,  showing  that  a  DC  current  is  generated.  The  induced  DC  current  will  eventually  heat  up  the 
diamond  lattice  due  to  electron-phonon  interactions  [35]  on  a  longer  time  scale,  which  we  neglected  in  the 
present  simulation. 

The  density  of  conduction  band  states  at  the  end  of  the  pulse  is  shown  in  Fig.4(a-c)  for  parallel  polar¬ 
ization  of  the  laser.  For  /  ~  1  TW/cm2  shown  in  Fig.4(a),  the  distribution  displays  isolated  Lorentzian 
peak  located  at  1.5  eV  above  the  conduction  band  minimum,  reflecting  the  fact  that  in  this  regime  transi¬ 
tions  occur  near  the  BZ  centre.  For  the  increased  field  intensity  in  Fig.4(b)  I  =  10  TW/cm2,  the  position 
of  the  low-energy  peak  is  unaffected,  however  the  Lorentzian  profile  is  distorted  by  the  laser  pulse  and  is 
supplemented  by  additional  satellite  structures  on  the  low  and  high  energy  sides  of  the  central  peak.  For 
the  higher  pulse  intensity  I  =  50  TW/cm2  shown  in  Fig.4(c),  the  density  of  conduction  band  states  extends 
due  to  release  of  electrons  with  higher  kinetic  energies.  Since  the  excess  kinetic  energy  of  the  photoexcited 
carriers  is  above  threshold  for  impact  ionization,  the  subsequent  energy  loss  during  binary  encounters  would 
limit  this  energy  distribution  within  the  bulk.  However,  since  large  fraction  of  electrons  have  already  been 
photoexcited  close  to  the  conduction  band  minimum,  the  scattering  probabiltiy  into  these  low-energy  states 
is  reduced  due  to  Pauli  blocking.  Thus  for  high  level  of  electronic  excitation  the  probability  for  impact 
ionization  should  become  unlikely,  such  that  electrons  can  gain  sufficiently  high  energies  from  the  laser 
field. 

In  Fig. 5  we  show  the  photoionization  yield  as  a  function  of  the  pulse  intensity  I  for  parallel  polarization 
of  the  electric  field.  In  the  low  intensity  regime  I  <  1  TW/cm2,  the  electron  yield  increases  linearly 
with  the  increase  of  I  due  to  highly  suppressed  one -photon  ionization  (with  negligible  final  densities  ne  < 
10  l4cr/F  :>>).  The  photoionization  probability  increases  steeply  with  the  increase  of  the  intensity  in  the  range 
1=1  TW/cm2  to  I  =  2  TW/cm2. 
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The  slope  of  the  theoretical  curve  can  be  fitted  to  n  ~  IN  scaling  law  over  that  range,  with  N  =  9 
corresponding  to  9-photon  ionization,  showing  the  deviation  from  the  perturbative  scaling  law  I5  occurs, 
pointing  out  that  electron  promotion  across  the  diamond  bandgaps  is  sensitive  to  the  detailed  characteristics 
of  the  substrate  band-structure  (cf.  also  [37,  38]).  For  I  «  3  TW/cnr  ionization  becomes  unlikely,  as  the 
electron  yield  decreases  by  one-order  of  magnitude.  That  is  because  interband  transitions  occur  away  from 
the  T  point,  where  the  wide  bandgaps  of  diamond  (~  10  eV)  suppress  ionization.  For  increased  light  pulse 
intensities  I  >  5  TW/cm2,  the  slope  of  the  curve  flattens  out  as  function  of  I.  In  this  high- intensity  regime, 
the  momentum  distribution  of  ionized  electrons  becomes  more  diffuse  and  extends  over  the  entire  Brillouin 
zone.  When  I  =  50  TW/cm2  the  earner  densities  are  ne  >  1020cm-3  showing  that  dielectric  breakdown 
occurs  [39].  Indeed  for  such  high  conduction  electron  density,  the  corresponding  bulk  plasma  frequency 
Up  =  (4 irrie)1/2  matches  the  laser  frequency  u>l,  such  that  electrons  can  absorb  very  efficiently  energy 
from  the  light  pulse  via  excitation  of  bulk  plasmon  resonance. 

The  non-perturbative  I9  scaling  law  in  the  low-intensity  regime  I  ~  1  TW/cm2,  can  be  expalined  by  we 
considering  a  simplified  two-band  model  taking  into  account  the  three-fold  degeneracy  of  the  valence  band 
at  the  T  point,  i.e. 

idtaVi  =  A (t)  ■  pViCac(t ) 

3 

idtac  =  Aac  +  E  A«  -  P evi  O-Vi  (/) 

2—1 

and  A  =  7  eV  is  the  bandgap.  Direct  interactions  among  valence  bands  vanish  because  of  the  even  sym¬ 
metry  of  their  wave-functions  under  spatial  inversion.  The  equation  of  motion  for  the  photoionization 
amplitude  can  be  integrated  formally  with  the  result 

ac(t)  =  -i  ^  f  dt' A{t')  ■  pCViaVi{t')  (13) 

i= l 

such  that  the  time-evolution  of  the  valence  band  amplitudes  becomes 


showing  that  the  Stark  mixing  has  introduced  indirect  dipole-dipole  couplings  among  valence  bands  for 
v.j  /  v j.  In  Fig.  6(a)  we  plot  the  time-evolution  of  the  absolute  value  of  the  Fourier  amplitudes  (for  laser 
intensity  I  =  1.5  TW/cm2)  by  neglecting  the  effect  of  induced  dipole  couplings.  For  this  low  intensity 
of  the  laser  field,  the  valence  band  is  only  weakly  distorted  hence  photoexictation  into  the  conduction  is 
highly  unlikely.  In  contrast,  the  strong  indirect  interaction  among  Stark-mixed  valence  bands  enables  the 
photoionization  of  bulk  diamond  as  shown  in  Fig.  6(b). 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


11 


The  enhancement  of  the  photoionization  yield  is  also  exhibited  in  the  time-evolution  of  the  optically 
induced  interband  polarization  in  the  wake  of  the  pulse  in  Fig.  6(c).  Due  to  the  induced  dipole-dipole  inter¬ 
actions  between  valence  bands,  the  polarization  displays  rapid  time  oscillations  with  period  corresponding 
to  the  bandgap  energy  7  eV,  otherwise  it  follows  adiabatically  the  laser  electric  field.  Thus  by  neglecting  the 
induced  dipole  couplings,  the  Fourier  transform  of  the  interband  polarization  in  Fig.  6(d)  exhibits  character¬ 
istic  perturbative  behavior,  it  includes  1,3  and  5th  harmonics  with  rapidly  decreasing  intensities.  This  is  in 
accord  with  the  conventional  Keldysh  theory,  which  neglects  energy-level  degeneracies  in  the  substrate  band 
structure.  Noticeably  however,  laser-induced  interactions  among  valence  bands  produces  non-perturbative 
harmonics  in  the  Fourier  transform  of  the  induced  polarization  with  enhanced  5th  harmonic  relative  to  the 
3rd  one.  This  effect  is  substantial,  since  the  5th  harmonic  is  enhanced  by  nearly  2  orders  of  magnitude 
relative  to  the  conventional  perturbative  multiphoton  result.  Thus  dipole  couplings  between  degenerate  va¬ 
lence  bands  create  a  highly  polarizable  intermediate  state  that  ionizes  very  efficiently  under  the  light  pulse 
to  produce  the  free  electrons  in  the  conduction  band.  It  is  worth  to  note  that  similar  results  are  found  in  Ref. 
[40],  showing  that  a  strong,  few-cycle  optical  field  is  capable  of  transforming  silica  into  a  state  of  highly 
increased  polarizability. 

We  have  compared  the  theoretical  calculations  of  the  ionization  yields  ne  reported  in  Fig.  5  as  a  function 
of  intensity  I  and  the  measured  generated  charge  Q  as  a  function  of  the  energy  per  pulse  E  taking  into 
consideration  the  space-dependence  of  the  laser  intensity  in  proximity  of  the  beam  waist.  Since  the  charge 
diffusion  length  during  the  pulse  is  much  smaller  than  both  of  the  waist  diameter  and  of  the  wavepacket 
length,  we  can  consider  the  generated  charge  density  as  only  dependent  on  the  local  energy  density,  defined 
as  (j)(r,  z)  =  I (r,  z)tl,  with  r  the  distance  from  the  beam  axis  and  z  the  position  along  the  beam: 

/oo  poo 

dz  /  2nrdr  ■  ne  ( <j>(r ,  z))  (15) 

-oo  J  0 

For  a  gaussian  beam  profile  attenuated  by  a  factor  exponentially  decreasing  with  the  absoiption  length: 

x 

(16) 

where  Ao  is  the  wavelength  in  vacuum.  The  constant  wq  and  zr,  have  been  evaluated  measuring  the  size  of 
the  graphitized  region  of  a  polycrystalline  diamond  at  the  bulk  graphitization  threshold,  which  are  respec¬ 
tively  about  2-3  and  35  pm  .  The  extinction  coefficient  k  has  been  assumed  to  depend  on  the  plasma  density 
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ne,  along  with  the  refractive  index  n,  according  to  the  formulas: 

n2  +  k2  =  ng  —  (cjp/wi)2 

2  nn  =  u2/u >2L-^g,  (17) 

where  coPtL  are  respectively  the  plasma  and  the  laser  frequencies,  p  and  m*  are  the  mobility  and  the  effective 
mass  of  the  charge  carriers,  respectively,  and  g  ~  1  is  a  constant,  depending  on  the  scattering  mechanism 
[41].  In  Eq.  17  a  great  importance  is  assumed  by  the  mobility  factor  p  decreasing  by  two  orders  of 
magnitude  [42]  when  ne  increases  from  1016  to  1019  cm  3.  As  a  consequence,  when  the  generated  charge 
density  reaches  the  order  of  1019  cm-3,  the  absorption  length  becomes  smaller  (possibly  much  smaller) 
than  zr,  and  the  laser  pulse  is  absorbed  long  before  reaching  the  region  of  the  lase  waist.  At  these  densities, 
the  plasma  frequency  is  as  low  as  between  1/7  to  1/3  of  ujr,  and  the  refractive  index  differs  negligibly  from 
its  static  value  2.42.  This  justifies  the  assumption  of  negligible  reflection-diffraction  effects  implied  by 
Eq.  16.  Note  that  ne  depends  on  <i>  in  Eq.  15,  and  6  depends  on  the  density  ne  via  the  extinction  coefficient 
in  Eq.  16.  Self-consistency  is  ensured  if,  along  each  ray  path,  we  calculate  values  of  0(r,  z)  for  increasing 
values  of  taking  into  account  the  previously  calculated  values  of  n(r' ,  z')  (and  thus  ne,  and  <f>)  for  z'  <  z. 

Fig.  7  shows  the  final  comparison  of  theory  with  experimental  data.  The  three  curves  reported  corre¬ 
spond  to  the  charge  calculated  from  Eq.  15  with  wq=2,  2.5  and  3  pm,  zr- 35  /ini.  The  Q(0)  values  derived 
from  experiment,  are  also  shown  together  with  the  theoretical  charge  values  calculated  assuming  a  carrier 
density  equal  to  the  ionization  yield  ne  in  Fig.  5,  uniformly  distributed  in  an  ellipsoidal  volume  element  of 
axis  wq  =  2.5  /mi  and  zr  =  35  pm. 


D.  Conclusion 

We  investigated  theoretically  the  photo-ionization  of  monocrystalline  CVD  diamond  subject  to  intense 
femtosecond  laser  irradiation  and  compared  the  simulations  with  experimental  measurements  that  were 
carried  out  under  the  very  same  conditions  used  for  preparing  three-dimensional  devices  by  fabricating 
conductive  columnar  electrodes  via  laser  bulk  graphitization.  The  theoretical  model  can  explain  the  exper¬ 
imental  data  starting  from  low  pulse  energies  and  reaching  about  E=  0.2  /J,with  a  focus  area  of  about  10 
pm2.  For  higher  intensities,  up  to  E= 5  pJ,  the  theory  correctly  explains  the  experiment  when  pulse  absorp¬ 
tion  due  to  the  high  plasma  density  is  phenomenologically  taken  into  account.  The  agreement  of  the  data 
with  the  theory  extends  over  six  orders  of  magnitude  in  terms  of  generated  charge  and  three  orders  of  mag¬ 
nitude  in  terms  of  laser  intensity.  The  graphitization  threshold  is  not  reached  in  the  bulk  because  non-linear 
absorption  limits  the  intensity  of  the  laser  field  and  energy  is  deposited  before  reaching  (i.e.  upstream  of) 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


13 


the  focus.  On  the  contrary,  when  we  prepare  the  devices  the  laser  is  initially  not  focused  inside  the  diamond 
bulk  but  the  conductive  column  is  started  from  a  localized  region  in  the  diamond  surface  where  surface 
defects  act  as  seeds.  In  order  to  fully  explain  the  graphitization  process  it  is  mandatory  to  include  a  model 
for  the  defects  involved.  Our  theoretical  approach  neglects  electron-electron  collisions  on  a  femtosecond 
timescale  which  eventually  cause  the  electron  avalanching  in  diamond.  Although  all  these  complex  mecha¬ 
nisms  are  not  fully  included  in  the  simulations,  we  find  that  the  experimental  data  is  quantitatively  described 
within  our  model  incorporating  realistic  bulk  band  structure  of  diamond.  Moreover  we  find  that  the  Stark 
effect  creates  a  highly  polarizable  transient  state  of  bulk  diamond  that  ionizes  very  efficiently  under  the  light 
pulse.  A  complete  understanding  of  the  physics  involved  is  required  to  tune  the  experimental  parameters 
for  preparing  devices  that  can  be  used  for  particle  detection  and  dosimetry. 

HI.  ULTRAFAST  PHOTOIONIZATION  AND  EXCITATION  OF  SURFACE-PLASMON-POLARITONS 

ON  DIAMOND  SURFACES 

A.  Introduction 

Experimental  ultrafast  SPP  photoexcitation  on  diamond  surfaces  was  not  achieved  yet,  even  though  their 
potential  imprinting  in  surface  relief  in  the  form  of  polarization-dependent  laser-induced  periodical  surface 
structures  (LIPSS,  surface  ripples)  was  numerously  evidenced  [43-45].  Such  experimental  studies  were 
devoted  to  the  design  and  fabrication  of  bio-sensors,  employing  the  biocompatibility  of  the  material,  by 
ablative  surface  nanostructuring  of  its  surface  with  high-intensity  femtosecond  (fs)  laser  pulses,  assuring 
precise  delivery  of  energy,  while  precluding  collateral  thermal  effects. 

Generally,  spatial  LIPSS  periods  A  are  known  to  depend  on  the  laser  wavelength  A  and  the  polarization 
of  the  laser  electric  field  e  and  the  number  of  laser  pulses  [46-50,  52-54],  The  surface  ripple  period  can  be 
slightly  less  than  A,  resulting  from  the  in-plane  weak  interference  of  the  incident  transverse  fs-laser  wave 
and  nearly  transverse  surface  polaritons  [55].  These  surface  electromagnetic  modes  are  photoexcited  by  the 
fs-laser  pump  pulse  via  its  scattering  on  permanent  or  laser- induced  (e.g.,  phase  transition  from  diamond 
to  glassy  or  diamond-like  carbon  phase)  cumulative  surface  relief  roughness  [50,  52-54],  or  prompt  laser- 
induced  ’’optical  roughness”  [56],  if  the  condition  5Pe[em]  <C  — 9?e[ed]  is  fulfilled  [49,  57]  for  the  dielec¬ 
tric  permittivities  of  the  photoexcited  surface  and  the  unperturbed  one.  .Meanwhile,  in  the  corresponding 
spectrally-naiTow  surface  plasmon  resonance,  occurring  for  the  photoexcited  surface  at  ’Re[em]  =  —  3?e[erf], 
the  short-wavelength,  longitudinal  surface  plasmons  can  similarly  interfere  with  the  incident  wave  or  among 
themselves  (for  counter-propagating  quasi-monochromatic  surface  plasmons),  inducing  surface  ripples  with 
periods  much  lower  than  A  (  A/2,  A/6,  ..)  [49,  58-60]. 
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Since  the  prompt  dielectric  permittivity  of  the  photoexcited  surface  appears  to  be  crucial  for  excitation 
of  either  near-wavelength  surface  polaritons,  or  sub-wavelength  surface  plasmons,  prompt  photoionization 
of  diamond,  directly  affecting  its  dielectric  permittivity,  should  be  explored  in  details.  There  are  numerous 
semi-empirical  approaches  to  explain  LIPSS  formation  e.g.  [47,  61-64],  corroborating  the  experimental 
evidence,  but  no  genuine  microscopic  approach  is  invoked  so  far.  The  basic  physical  processes  involve  ex¬ 
citation  of  electron-hole  pairs,  often  parameterized  by  Keldysh  approximate  formulas.  Photoionization  may 
produce  highly  energetic  electrons  that  collisionally  ionize  the  valence  band  and  produce  more  electrons  in 
the  conduction  band.  The  multiplication  of  carriers  may  cause  optical  breakdown  of  bulk  diamond.  The  col¬ 
lective  response  of  charge  carriers  screens  out  the  laser  electric  field  inside  the  bulk  when  the  number  density 
is  sufficiently  large.  At  some  instant  of  time  the  bulk  dielectric  function  may  become  negative  at  the  laser 
wavelength,  allowing  excitation  of  SPP  at  the  rough  surface  and  LIPSS  formation  via  the  optical  interference 
mechanism.  The  dielectric  properties  of  the  laser-irradiated  material  in  most  cases  are  parameterized  with 
Drude  model  [52,  59,  65-68],  which  combines  the  ground  state  response  with  the  laser-induced  free-carrier 
response.  This  model  usually  requires  three  free  parameters  -  the  number  density  of  electron-hole  pairs,  the 
free-carrier  effective  mass  and  the  Drude  damping  time,  which  are  adjusted  to  fit  experimental  data.  Ref. 
[69,  70]  proposed  more  elaborate  model  for  the  optical  dielectric  function,  which  implements  state-  and 
band-filling  effects,  renormalization  of  the  band  structure  and  free-carrier  response.  The  dielectric  function 
of  laser-excited  silicon  was  studied  from  first  principles  using  the  time-dependent  density  functional  theory 
(TDDFT)  [71].  A  distinguishing  feature  in  the  linear  response  of  the  photoexcited  silicon  is  a  plasmon  peak 
with  large  Drude  damping  time  as  short  as  re  ~  1  fs,  despite  the  fact  that  collisional  effects  are  neglected  in 
the  TDDFT  simulation.  The  real  part  of  dielectric  function  was  well  fitted  by  a  Drude  free-carrier  response 
showing  that  Re  [cm]  is  sensitive  to  the  total  number  density  of  excited  electrons  and  not  to  the  detailed 
distribution  of  electron-hole  pairs,  while  sensitivity  to  the  nonequilibrium  distribution  of  the  phototexcited 
carriers  manifests  in  the  imaginary  part  of  the  dielectric  constant.  Subsequently,  TDDFT  was  applied  to 
study  ablation  of  silica  subjected  to  ultrashort  laser  pulses  [72],  The  comparison  between  the  estimated  sur¬ 
face  ablation  threshold  and  the  experimental  data  suggests  a  non-thermal  mechanism  in  the  laser  ablation  of 
silica  by  fs-laser  pulses,  furthermore  theoretical  ablative  crater  depths  agree  with  the  measured  ones.  The 
drawback  of  this  approach  is  its  limitation  to  very  short  laser-matter  interaction  timescales  (less  than  10  fs). 

We  present  theoretical  results  for  the  laser  ablation  and  subsequent  LIPSS  formation  on  diamond  sur¬ 
faces  subjected  to  normally  incident  515-nm,  200-fs  laser  pulses  due  to  the  excitation  of  SPPs.  Our  the¬ 
oretical  modeling  of  LIPSS  formation  on  diamond  surfaces  is  based  on  numerical  solution  of  the  time- 
dependent  Schrodinger  equation  (TDSE)  in  bulk  diamond  subjected  to  a  single  intense  laser  pulse.  The 
theory  describes  the  electron  dynamics  quantum  mechanically  in  the  single-active-electron  approximation. 
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Collisional  de-excitation  of  the  photoexcited  carriers  and  subsequent  impact  ionization  are  treated  within 
rate  equation  approach  and  an  optical  breakdown  threshold  is  derived.  Due  to  the  contribution  of  the  impact 
ionization  the  real  part  of  the  bulk  dielectric  constant  of  the  irradiated  diamond  becomes  negative  in  the 
trailing  edge  of  the  pulse  resulting  in  plasma  that  is  opaque  to  the  incident  radiation.  The  inhomogeneous 
energy  deposition  in  the  surface  was  modeled  with  the  Sipe-Drude  efficacy  factor  theory  [64,  67]  in  terms  of 
time-dependent  dielectric  function  of  free  carriers.  The  applicability  of  this  efficacy  factor  theory  for  LIPSS 
formation  in  laser-irradiated  dielectrics  was  confirmed  by  numerical  solutions  of  the  Maxwell’s  equations 
at  statistically  rough  surfaces  [73]. 


B.  Methods,  Assumptions  and  Procedures 


In  order  to  model  theoretically  LIPSS  formation  in  femtosecond-laser-excited  diamond,  we  apply  the 
ah  initio  theory  developed  by  Sipe  [47].  In  this  picture,  the  laser  beam  is  incident  on  a  rough  surface,  the 
(permanent  or  laser-induced)  roughness  is  assumed  to  be  confined  within  a  surface  region  (selvedge)  of 
thickness  l  much  smaller  than  the  laser  wavelength  A.  The  optically-induced  polarization  in  the  selvedge 
generates  surface-scattered  waves  that  interfere  with  the  refracted  laser  beam  leading  to  inhomogeneous 
energy  deposition  into  the  surface.  The  inhomogeneous  energy  absorption  can  be  described  by  the  function 


A(k)  ~  |&(k)|7?(k;  Ki), 


(18) 


where  Kly  is  the  component  of  the  laser  propagation  wave  vector  parallel  to  the  surface,  b(n)  is  a  measure  of 
surface  roughness  at  wave-vector  k  =  ( kx ,  ny)  and  //(k,:  k,  )  is  an  efficacy  factor  describing  the  contribution 
to  the  energy  absoiption  at  the  LIPSS  wave  vector  k.  The  prediction  of  Eq.18  is  valid  if  the  selvedge  thick¬ 
ness  is  small  compared  to  the  LIPSS  period,  i.e.  kI<1  should  be  satisfied.  The  efficacy  factor  essentially 
incorporates  the  modification  of  the  surface  morphology  and  the  variation  of  the  dielectric  constant  e  of  the 
photoexcited  diamond.  For  normally  incident  s-polarized  laser  pulse  with  k,  =  0,  the  efficacy  factor  (as  a 
function  of  the  normalized  wave-number  k  =  A/A)  can  be  written  as  t/(k)  =  47r|9te[p(#e)]|  with 


o(k)  = 


(19) 


where  the  response  functions  hss  and  hKK 


M«)  =  2i-^,  /.„(*)  =  2  ,  (20) 

hvy  i  Him  EKy  Him 

are  given  in  terms  of  the  transient  bulk  dielectric  function  e(uj;t)  ,  kv  =  V n2  —  1  and  nm  =  V k2  —  e, 
the  Fresnel  transmission  coefficient  ts  =  2/(1  +  \Je  —  1)  in  the  absence  of  the  selvedge,  the  effec¬ 
tive  transverse  susceptibility  function  =  (e  —  l)/47r{e  —  (1  —  /)(e  —  l)[/i(s)  —  Rhi(s)]},  the  surface 
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roughness  characterized  by  shape  s  and  filling  /  factors,  R  =  (e  —  l)/(e  +  1)  and  the  shape  functions 
h(s)  =  V s2  +  1  —  s,hi(s )  =  (Vs2  +  4  +  s)/2  —  Vs2  +  1.  When  3?e[e]  <  0,  the  response  function  /iss 
exhibits  small  kinks  near  the  light  line  k  ~  1,  in  contrast  hKK  exhibits  sharp  resonance  structure  due  to  the 
excitation  of  surface  plasmons  and  diverges  at  the  (complex)  SPP  wave-number  kspp  =  u)/cy/e/(\  +  e). 

As  in  section  II  photoexcitation  and  the  dielectric  response  of  laser-irradiated  diamond  are  treated  in 
independent  particle  approximation  based  on  the  3D  TDSE.  In  a  long-wavelength  approximation  the  light 
pulse  is  represented  by  a  spatially  uniform  time-dependent  electric  field  and  velocity  gauge  is  used  through¬ 
out  the  calculations  [74].  The  static  bulk  band  structure  is  represented  by  the  lowest  4  valence  bands  and 
16  unoccupied  conduction  bands.  The  Brillouin  zone  was  sampled  by  a  Monte  Carlo  method  using  2000 
randomly  generated  k-points.  The  time  step  for  integration  of  the  equations  of  motion  was  St  =  0.03  a.u. 

For  the  static  band  structure  shown  in  Fig.l  earner  excitation  occurs  through  the  direct  gap  at  the  T 
point,  however  excitation  into  higher  lying  conduction  bands  is  also  a  relevant  process  for  the  considered 
laser  intensity  range  lG|  1 ,50]  TW/cm2. 

During  the  irradiation  of  the  diamond  surface  with  pulsed  200  fs-laser,  the  total  number  of  electrons 
generated  into  the  conduction  band  is  given  by  a  Brillouin  zone  integral 


Pe(t)  =  f 

VVM 4?r 


(21) 


where  /nk  is  the  occupation  number  of  the  n-th  conduction  band  and  k  is  the  crystal  momentum.  The 
electronic  excitation  energy  per  unit  cell  is  given  by 


Ef 


■H)  =  £  / 

.„«A 


d3k 

-j— 3  (V’nk(f)l^tlV’nk(f))  -  Eo, 
bz  4vr- 


(22) 


where  V'nk(f)  are  the  time-evolved  Bloch  orbitals  of  valence  electrons  and  Eq  is  the  ground-state  energy. 
The  time  evolution  of  the  free-electron  density  is  shown  in  Fig.  8a,  for  linearly  polarized  electric  field  along 
the  (1,1,1)  direction  with  intensity  30  TW/cm2.  Carrier  generation  occurs  efficiently  prior  to  the  peak  of 
the  pulse.  Transient  charge  density  oscillations  following  the  laser  period  are  due  to  quiver  motion  of  free 
electrons  in  the  electric  field.  An  electron-hole  plasma  (EHP)  with  number  density  exceeding  102 1  cm~ 3 
is  established  shortly  after  the  peak  intensity.  The  cycle-averaged  photoelectron  yield,  shown  in  Fig.  8b, 
is  a  slowly  varying  function  of  time.  Carrier  generation  on  the  rising  edge  of  the  pulse  competes  with 
recombination  on  the  trailing  edge  of  the  pulse  to  determine  the  final  photoionization  yield.  Recombination 
of  carriers  becomes  unlikely  with  the  increased  laser  intensity,  cf.  also  Fig.  8b.  The  cycle  averaged  electron 
yield  includes  contributions  due  to  creation  of  real  as  well  as  virtual  electron-hole  pairs.  Since  adiabatic 
evolution  does  not  produce  any  real  excitation  of  the  crystal,  the  carrier  density  should  be  calculated  with 
respect  to  adiabatically  evolved  ground  state  orbitals  that  are  obtained  from  the  static  Bloch  orbitals  with 
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shifted  crystal  momentum  k(f)  =  k  +  A (t),  i.e.  p(t)  =  pe(t)  —  pad(t),  where  the  adiabatic  density  is 

f  d3  k  o 

Pad{t )  =  ^  ^  /  T  3  K^nkl^nkp))  I  j  (23) 

£n<0,/BZ  7F 

and  {|?Tnk}}  are  the  lattice-periodic  Bloch  states.  The  number  density  of  photoexcited  earners  is  shown 
in  Fig.  8c  .  It  can  be  seen  that  discarding  contributions  of  virtually  excited  electron-hole  pairs  leads  to 
reduction  in  the  number  density  by  an  order  magnitude  near  the  peak  of  the  pulse.  The  electronic  excitation 
energy  is  shown  in  Fig.  9(a)  for  laser  intensity  30  TW/cnr.  The  temporal  variation  of  the  cycle-averaged 
energy  gain  follows  closely  the  envelope  of  the  laser  pulse  during  the  first  half  of  the  driving  pulse  and 
reaches  1.5  eV/atom  at  the  peak  of  the  pulse  that  is  small  as  compared  to  the  cohesive  energy  of  diamond 
7.37  eV/atom.  After  the  pulse  peak,  electron-hole  pairs  recombine  by  transferring  part  of  their  energy  back 
to  the  radiation  field.  Energy  exchange  is  not  completely  reversible  since  the  time  delay  in  restoration 
of  equilibrium  gives  rise  to  a  net  energy  gain  of  0.5  eV  per  carbon  atom  after  the  end  of  the  pulse.  The 
deposited  energy  increases  steadily  with  the  increase  of  the  intensity,  i.e.  for  1=50  TW/cm2,  it  reaches  2 
eV/atom.  Since  this  excitation  energy  is  still  lower  than  the  diamond  cohesive  energy,  Fig.  9(b)  shows  that 
ablation  threshold  is  not  reached  up  to  1=50  TW/cm2. 

For  the  200fs  pulse  duration  and  intensities  lower  than  50  TW/cm2  the  electron  density  produced  by 
photoionization  is  below  the  critical  one.  That  suggests  that  impact  ionization  is  the  relevant  process  that 
determines  the  optical  breakdown  threshold.  In  Fig.  10  (a-c)  we  plot  the  density  of  conduction  states  after 
the  end  of  the  pulse.  It  is  seen  that  the  laser  has  created  electron-hole  pairs  with  well-defined  energies.  This 
non-thermal  distribution  relaxes  towards  the  equilibrium  Fermi-Dirac  distribution  on  a  time  scale  ranging 
over  few  tens  of  a  femtoseconds  to  a  picosecond  [75,  76]  without  changing  the  electron  number  density. 
Photoelectrons  are  excited  into  the  lowest  conduction  band  across  the  direct  gap  (with  energies  2  eV  above 
the  conduction  band  minimum)  and  substantial  fraction  of  carriers  occupy  higher  lying  conduction  bands 
with  energy  above  threshold  for  impact  ionization  (specified  by  the  indirect  gap  5.4  eV).  These  highly 
energetic  electrons  may  collisionally  de-excite  to  lower  energy  states  and  their  excess  energy  is  spent  to 
promote  valence  electrons  into  the  conduction  band.  Assuming  that  the  time  evolution  of  the  electron 
density  is  governed  by  a  rate  equation  [77,  78] 

^  =  G(t )  -  R(t)  +  wimp(I)p  (24) 


including  carrier  generation  G(t)  and  recombination  It  it)  rates  supplemented  by  an  intensity-dependent 
impact  ionization  rate  Wimp  obtained  as  a  weighted-average  of  the  field-free  ionization  rate 


/e°°  dep(e;  I)Pimp(e) 
fo°°  dep(e;I ) 


(25) 
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here  p(e;  I)  is  the  density  of  conduction  states  after  the  end  of  the  pulse  (cf.  Fig.  lOa-c), 

p(e'i1)=^2  /  X-3/nk(-f)<5(e-enk),  (26) 

en>OjBZ  n 

Pirnp(f)  =  -Fb(e  —  e,;)4"5  is  the  energy-dependent  impact  ionization  rate  for  diamond,  et  is  the  threshold  for 
impact  ionization  (5.42  eV)  and  Pq  =  3.8  x  1010s~1el4_4'5  [79]. 

In  contrast  to  the  standard  perturbative  result  based  on  Keldysh  theory  valid  for  monochromatic  laser 
radiation  the  calculated  carrier  generation  and  recombination  rates  shown  in  Fig.  11a  do  not  follow  the 
temporal  profile  of  the  laser  pulse.  This  result  suggests  that  the  pulse  shape  and  pulse  duration  are  relevant 
control  parameters  for  non-adiabatic  electron  dynamics  in  the  laser  irradiated  diamond.  The  key  features 
are  generation  of  dense  plasma  50fs  prior  to  the  pulse  peak  and  subsequent  laser-induced  recombination  of 
electron-hole  pairs  in  the  trailing  edge  of  the  pulse. 

Fig. lib  shows  the  impact  ionization  rate  that  depends  in  highly  non-linear  way  on  the  laser  intensity. 
This  non-linear  and  non-monotonic  intensity-dependence  reflects  the  population  of  higher-lying  conduction 
bands  (cf.  also  Fig.  10).  It  is  seen  that  the  impact  ionization  rate  reaches  few  tens  of  inverse  picosecond  for 
I  >  15  TW/cnr.  In  Fig.  12  we  plot  the  EHP  density  with  and  without  the  impact  ionization  term.  This 
comparison  demonstrates  that  photoionization  produces  the  seed  electrons  needed  for  the  impact  ionization 
on  the  rising  edge  of  the  pulse  and  then  the  conduction  electron  density  grows  exponentially  after  the  pulse 
peak  resulting  in  dense  plasma  (with  density  1022  cm~3)  50fs  after  the  peak  of  the  pulse. 

Since  the  absorption  of  the  femtosecond  laser  pulses  in  diamond  results  in  the  generation  of  nearly  free 
electrons  in  the  conduction  band  on  timescales  smaller  than  the  electron-phonon  relaxation  time  [80],  we 
describe  the  linear  response  of  the  photoexcited  diamond  by  a  free-carrier  Drude  response  [67]  using  of 
time-dependent  plasmon-pole-approximation  for  the  density-density  correlation  function  of  the  Coulombi- 
cally  interacting  electron  gas  [81] 

S(t,  t')  =  —  9(t  —  t')w3/2(f')w”1//2(t)  sin  f  drLOp(r),  (27) 

J? 

where  6(t)  is  the  Heaviside  step  function,  up(t)  =  (  j  is  the  bulk  plasma  frequency,  Eq  and  me  are 
the  vacuum  permittivity  constant  and  the  free-electron  mass,  respectively.  In  long  wavelength  approxima¬ 
tion  the  spatial  dispersion  of  the  bulk  plasmon  is  neglected.  The  Fourier  transformation  of  the  correlation 
function  is  the  transient  frequency  dependent  inverse  dielectric  function  of  the  free-electron  plasma 

t 

£-\uj-t)  =  1+  J  (28) 

—  OO 

where  6  =  l/re  is  a  free-carrier  polarization  dephasing  rate,  which  we  shall  treat  as  a  free  parameter.  If  the 
time  delay  in  the  build  up  of  screening  in  the  optically  excited  plasma  can  be  neglected,  the  classical  Drude 
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dielectric  function  is  recovered  e~ 1  (uj:  t )  =  a?2/ (a;2  —  cc2(t))  with  parametric  time  dependence  of  the  bulk 
plasma  frequency. 

In  Fig.  1 3a-c,  we  plot  the  real  and  imaginary  parts  of  the  dielectric  function  for  laser  intensity  1 8  TW/cnr . 
The  screening  charge  density  accumulates  during  the  first  half  of  the  pulse  (t  <  0).  Over  that  time  interval 
the  laser  frequency  is  above  the  plasma  frequency  and  the  diamond  surface  remains  transparent  to  the 
incident  radiation.  The  frequency  dependent  dielectric  function  displays  oscillations  in  the  spectral  range 
below  the  laser  frequency  due  to  the  time  lag  in  the  build  up  of  screening.  Because  of  the  impact  ionization, 
the  laser  frequency  falls  off  below  uip  after  the  pulse  peak  (t  =  25  fs)  when  the  plasma  is  reflective  for 
the  incident  radiation  and  an  optical  breakdown  threshold  is  reached.  In  this  regime,  the  dielectric  function 
essentially  exhibits  the  Drude  form  with  time-dependent  bulk  plasma  frequency  iop(t).  For  the  transiently 
increasing  carrier  density,  5Re[c]  passes  the  narrow  surface  plasmon  resonance  at  5Re[e]  =  —1,  with  kspp  A> 
1  and  A/A  <C  1,  and  becomes  large  and  negative  in  the  trailing  edge  of  the  pulse  (t  >  100  fs)  with 
Ae  [s]  >  0,  with  corresponding  nspp  >  1.  During  this  plasmonically-active  phase  of  the  laser-irradiated 
diamond  the  SPP-laser  interference  mechanism  of  inhomogeneous  energy  deposition  is  effective  and  leaves 
permanent  imprints  on  the  surface  morphology  after  the  conclusion  of  the  pulse. 


C.  Results  and  Discussion 

To  make  possible  identification  and  interpretation  of  experimentally  obtained  SPP  modes  we  plot  the 
efficacy  factor  as  a  function  of  the  wave  vector  k,  in  a  narrow  laser  intensity  range  above  the  optical  break¬ 
down  threshold  in  Fig.  14  a-b.  The  transient  bulk  dielectric  function  was  evaluated  at  the  laser  wavelength, 
i.e.  eit)  =  The  surface  roughness  was  modeled  as  a  collection  of  spherically-shaped  islands 

corresponding  to  standard  values  s  =  0.4  and  /  =  0.1  for  the  shape  and  tilling  factors  respectively.  For  nor¬ 
mally  incident  light  pulse,  the  numerical  results  are  weakly  dependent  on  the  specific  parameters  describing 
surface  morphology  and  therefore  the  transient  dielectric  constant  is  the  most  significant  in  determining  the 
efficacy  factor.  Here  we  demonstrate  that  the  main  features  in  the  inhomogeneous  energy  deposition  in  the 
surface  as  represented  by  the  position  of  the  peaks  of  the  transient  efficacy  factor  are  in  correspondence  with 
the  experimentally  observed  LIPSS  periods. 

In  a  very  narrow  laser  intensity  range,  when  the  laser  frequency  nearly  matches  the  surface  plasma 
frequency  (Fig.  14a),  the  efficacy  factor  has  large  contribution  due  to  excitation  of  the  surface  plasmon 
resonance  (SPR).  In  this  early  stage,  the  spatial  extension  of  the  electromagnetic  field  inside  the  bulk  asso¬ 
ciated  with  SPP  is  determined  by  the  SPR  decay  constant  Krn,  which  at  short  wavelengths  nm  — >  I  //  defines 
a  skin-depth  ls  =  1/ nm  ~  l  leading  to  strong  concentration  of  the  electromagnetic  field  in  the  thin  selvedge 
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region.  SPPs  need  finite  time  to  build  up  to  incorporate  the  details  of  the  surface  relief  and  interfere  with  the 
laser  to  modify  the  Fourier  components  of  the  surface  roughness  function  6(k)  via  periodic  laser  ablation. 
In  this  regime  the  deposition  of  laser  energy  into  the  surface  plasmon  wavevector  causes  formation  of  fine 
ripples  with  spatial  periods  around  100  nm,  as  observed  in  the  periphery  of  the  ablative  craters,  cf.  Fig.  16. 
The  transverse-magnetic  characteristic  of  the  SPP  determines  the  orientation  of  the  surface  ripples.  At  a 
later  time,  the  transiently  increasing  number  of  conduction  electrons  makes  the  dielectric  constant  large  and 
negative  at  the  laser  wavelength,  the  intensity  map  of  t](k:  0)  shrinks  and  concentrates  on  the  outer  part 
of  the  circle  k  =  1  (cf.  Fig  14b)  which  clearly  can  be  associated  with  the  formation  of  near- wavelength 
surface  ripples  oriented  perpendicularly  to  the  laser  polarization.  At  the  longer  wavelengths  with  k  — >•  0 
the  skin  depth  ls  =  c/lop  ~  80  nm  is  much  smaller  than  the  laser  wavelength.  Therefore,  above  the  SPR 
excitation  threshold,  the  transiently  increasing  carrier  density  results  in  a  shift  of  the  SPP  wave  number 
from  the  high  spatial  frequency  region  towards  the  light  line  (also  causing  expansion  of  the  skin  depth),  and 
this  red  shift  is  highly  sensitive  on  the  carrier  density  (or  laser  intensity),  cf.  Fig.  15a.  The  surface  plasmon 
peak  in  the  efficacy  factor  is  also  affected  by  the  relaxation  time  parameter  re  as  shown  in  Fig.  15b  from 
the  experiments  described  in  .  If  re  is  decreased  to  10  fs,  the  surface  plasmon  cusp  turns  into  a  dip,  which 
hinders  the  efficient  energy  absorption  at  the  surface  plasmon  wavevector.  This  dependence  suggests  that 
the  Drude  carrier  relaxation  time  parameter  influences  prompt  feedback  mechanisms  involved  in  the  forma¬ 
tion  of  surface  ripples.  Indeed  in  the  high-frequency  limit  with  LdfTe  ^$>  1,  the  metalized  surface  behaves 
as  a  nearly  ideal  inductor,  while  in  the  low-frequency  limit  c 0LTe  <C  1,  the  resistive  Ohmic  losses  result  in 
electron  heating  in  the  skin  layer. 

Because  the  efficacy  factor  theory  does  not  fully  account  for  interpulse  feedback  processes  that  are 
undoubtedly  important  in  the  detailed  development  of  morphological  features  on  the  diamond  surface  our 
theoretical  results  are  not  directly  applicable  to  the  multipulse  phase  of  LIPSS  formation.  Flowever  the 
experimental  data  shows  that  once  surface  ripples  are  formed,  exposure  by  subsequent  pulses  has  little 
effect  on  their  spatial  period  and  location,  thus  LIPSS  formation  should  be  possible  already  for  a  single¬ 
pulse  irradiation,  provided  that  SPP  can  be  excited,  e.g.,  by  surface  defects  [83].  Once  LIPSS  are  formed, 
the  spectrum  of  the  surface  roughness,  b(n)  contains  peaks  at  the  SPP  wavenumber  causing  enhanced 
inhomogeneous  energy  deposition  and  further  growth  of  LIPSS  as  is  also  evidenced  from  the  SEM  images  in 
Fig.l6a-b.  Furthermore  the  subsequent  pulses  interact  with  periodically  structured  surface  hence  a  grating- 
assisted  laser-surface  coupling  becomes  effective  [49]  causing  a  decrease  of  the  ripple  wavelength.  The 
experimental  data  shows  only  minor  modification  of  LIPSS  periods  that  is  consistent  with  the  hypothesis 
in  Ref.  [49]  that  because  of  strong  thermal  effect  at  the  crater  center,  the  grating-assisted  coupling  is  weak 
and  the  ripple  wavelength  is  unaffected  by  higher  exposure,  i.e.  depends  weakly  on  the  superimposed  pulse 
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number. 


D.  Conclusion 

In  this  section  the  transition  of  diamond  into  plasmonically-active  material  and  the  subsequent  formation 
of  ripples  on  its  surface  caused  by  SPP-laser  interference  were  investigated  via  first  principle  approaches. 
The  numerical  simulations  of  periodic  laser  energy  deposition  on  photo-excited  diamond  surface  revealed 
the  formation  of  fine  ripples  oriented  perpendicularly  to  the  laser  polarization  in  a  narrow  range  of  in¬ 
tensities  above  the  optical  breakdown  threshold.  The  transient  increase  of  the  carrier  density  above  this 
threshold  resulted  in  the  formation  of  near-wavelength  surface  ripples.  In  addition  we  demonstrate  that 
transient  changes  of  the  dielectric  permittivity  reflect  the  chronological  order  of  the  ripple  formation.  Semi- 
quantitative  agreement  with  the  experimentally  observed  periods  of  the  periodic  structures  was  obtained. 
The  realistic  insight  into  the  microscopic  processes  gained  by  our  combined  theoretical  and  experimental 
study  will  allow  to  predict  and  potentially  extend  the  capabilities  of  the  dielectric  nano-optical  meta  devices 
and  the  photo-modulated  metasurfaces  based  on  them. 


IV.  HIGH  HARMONIC  GENERATION  FROM  BULK  DIAMOND  DRIVEN  BY  INTENSE 

FEMTOSECOND  LASER  PULSE 

A.  Introduction 

High  harmonic  generation  is  a  secondary  process  of  laser  frequency  conversion  to  its  multiples  driven 
by  the  coupling  of  the  oscillating  electric  field  of  an  intense  ultrashort  laser  pulse  to  electrons  comprising 
matter.  The  microscopic  mechanisms  underlying  this  phenomenon  have  been  elucidated  experimentally 
demonstrating  that  high  intensity  laser  pulses  in  the  UV  and  IR  spectrum  generate  high  order  harmonics  in 
gaseous  and  molecular  targets  [8  -88],  plasmas  [89]  and  in  solid  materials  [16]  These  experiments  have 
shown  that  sources  based  on  high  harmonic  generation  (HHG)  have  spatial  and  temporal  coherence,  high 
directionality  and  polarization  but  due  to  the  low  efficiency  of  the  process,  the  intensities  of  the  generated 
harmonics  are  orders  of  magnitude  lower  than  that  of  the  exciting  field  hence  less  intense  than  EUV,  XUV 
and  soft  x-ray  radiation  produced  in  the  huge  facilities  such  as  FEL  accelerator  [90]  and  Synchrotron  Radi¬ 
ation  (SR)  [91].  Nevertheless  they  can  be  potentially  used  to  obtain  coherent  radiation  of  sub-femtosecond 
to  attosecond  duration  in  a  broad  wavelength  range  (XUV  to  X-ray)  [92-94]  capable  of  inducing  physical 
processes  on  very  short  time  and  space  scales  leading  to  applications  such  as  attosecond  electron  dynamics 
[95],  molecular  tomography  and  protein  crystallography  [96],  etc. 
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Experimental  studies  of  HHG  in  semiconductors  and  dielectrics  have  been  reported  recently,  e.g.  [16- 
18].  In  solid  crystals  the  interplay  between  generation  of  carriers  via  multiphoton  excitation  over  the 
bandgap  at  high  intensities  leading  to  catastrophic  breakdown  of  the  material  and  HHG  is  carefully  ex¬ 
plored  to  find  regimes  in  which  only  HHG  takes  place  avoiding  material  collapse.  In  the  pioneering  study 
of  Ghimire  et  al.  [  1 6]  long  wavelength  laser  radiation  is  used  to  suppress  multiphoton  ionization  and  to 
obtain  HHG.  Experimental  and  theoretical  study  of  HHG  spectra  from  rare-gas  solids  [17]  exposed  char¬ 
acteristic  features  of  HHG  in  solids:  abrupt  appearance  of  a  second  plateau  as  the  driving  laser  intensity 
increases  over  a  narrow  intensity  range  that  consists  of  a  few  harmonics;  complicated  dependence  of  the 
cutoff  of  both  plateaus  on  the  laser  field,  i.e.  following  neither  linear  (solid  phase)  nor  the  square -root  scal¬ 
ing  (gaseous  phase);  photon  energy  of  the  HHG  in  solids  exceeding  the  maximum  band  separation  between 
the  highest-valence  and  lowest-conduction  bands.  Studies  based  on  a  model  describing  electrons  tunneling 
from  the  valence  to  a  single  conduction  band  where  they  are  accelerated  by  the  driving  laser  field  and  experi¬ 
ence  Brgg  reflections  from  the  edge  of  the  Brillouin  zone  [97,  98]  Interband  transitions  as  well  as  interband 
coherences  between  conduction  bands  and  valence  bands  were  taken  into  consideration  in  other  models 
[99-101]  In  the  recent  theoretical  approaches  used  to  rationalize  HHG  in  crystalline  solids  including  HHG 
in  diamond  varying  from  Time  Dependent  Density  Functional  Theory  (TDDFT),  [19,  20,  102,  103]  Time 
dependent  Schrodinger  equation  (TDSE)  [100],  band  resolved  Semiconductor  Bloch  equations  [  01,  104] 
to  semiclassical  methods  [105]  the  latter  being  more  intuitive  than  the  former,  the  underlying  microscopic 
processes  are  still  explored  to  elucidate  further  the  dependence  of  HHG  in  solids  on  the  electronic  band 
structure  and  on  the  interplay  between  inter-  and  intraband  transitions  as  well  as  on  the  collective  behavior 
of  the  high  density  electrons  and  holes  coupled  to  the  lattice.  In  this  respect  methods  using  TDSE  including 
full  band  structure  treatment  together  with  saddle  point  methods  may  provide  very  consistent  and  thorough 
theoretical  treatment  of  the  HHG  generation  from  solid  phase.  In  view  of  the  fact  that  HHG  in  wide  band- 
gap  semiconductors  may  lead  to  brighter  sources  due  to  the  high  laser  induced  carrier  density  and  ability 
to  withstand  the  driving  laser  fields  necessary  for  inducing  strong  nonlinear  and  non  perturbative  effects 
leading  to  HHG  in  the  present  work  we  explore  the  feasibility  for  femto  and  attosecond  pulse  generation 
from  diamond  subjected  to  intense  15  fs  laser  pulse  with  mid-infrared  wavelength  by  exploring  further  and 
deeper  the  non-perturbative  scaling  of  the  generated  harmonics  with  the  laser  field  and  their  bandstructure 
dependence. 
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B.  Methods,  Assumptions  and  Procedures 


In  velocity-gauge,  the  TDSE  in  single-active  electron  approximations  is  given  by 

*^|^k(i))  =  (H0  +  Hiut(t))\i’v]i(t)}  (29) 

where  nk  labels  the  initially  occupied  valence  band  states  of  definite  crystal  momentum  k,  Hq  is  the  Hamil¬ 
tonian  of  the  undistorted  lattice 

H0  =  ±p2  +  V(r)f  (30) 

and 


Hhlt(t)  =  A (£)  ■  p  (31) 

is  the  interaction  Hamiltonian  between  the  laser  field  and  the  electron,  p  is  the  momentum  operator  and 
V (r)  is  the  periodic  ion-lattice  potential.  The  vector  potential  is  parametrized  by  a  Gaussian  function  as  in 
section  II.  The  eigenstates  of  the  field-free  Hamiltonian  arc  Bloch  states  of  definite  crystal  momentum 

H0\nk)  =  enk|nk).  (32) 


We  expand  the  time-dependent  states  over  this  basis 

1 1pm(t)}  =  ^2  lnk)(rak|  1pm(t)}  (33) 

nk 

and  calculate  the  optically-induced  electron  current  in  direction  of  the  laser  electric  field  is  given  by  the 
Brillouin  zone  integral 

m-L*™  (34) 

where 


<7(M)  =  X]^k(*)le  '  PlV’vkC*)), 


(35) 


the  diamagnetic  gauge  current  contribution  can  be  neglected,  because  it  does  not  contribute  to  generation  of 
high  harmonics.  The  HHG  spectrum  can  be  presented  as  a  coherent  sum 

d3  k 


Jbz  4vt3  Jk’“ 


I(u)  = 

over  contributions  from  recombination  of  electron-hole  pairs  with  wave-vector  k 

Jk,u>  =  [  dteiutJ(k,  t). 


(36) 


(37) 
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The  valence  electron  states  are  propagated  in  time  for  small  equidistant  steps  of  5  =  0.7  attoseconds  by 
applying  the  Cayley  transformation 

\ipv{t  +  S))  «  (l-i^H(t)\  +  i5-H(t)^j  \^v(t)}  (38) 

where  I  is  the  unit  matrix. 


C.  Results  and  Discussion 

Fig.  17(a-c)  shows  F1HG  spectra  for  three  laser  intensities  I  =  2, 10  and  30  TW/cm2.  For  the  low  inten¬ 
sity  shown  in  Fig.  17(a),  the  HHG  is  characterized  by  odd  harmonic  structure  due  to  the  inversion  symmetry 
of  the  diamond  lattice  and  consists  of  the  1st  and  3rd  sub-bandgap  harmonics  and  the  5th  harmonic.  Since 
the  direct  band  gap  at  the  T -point  is  nearly  7  e V,  the  lowest  excitation  of  valence  electrons  across  this  gap 
occurs  via  5  photon  absorption.  At  this  lower  intensity,  the  HHG  ends  abruptly  at  the  position  of  5th  har¬ 
monic,  such  that  the  high-energy  cutoff  is  determined  by  the  direct  gap.  For  the  increased  intensity  shown 
in  Fig.  17(b),  the  spectrum  extends  to  higher  harmonic  orders,  a  single  plateau  is  seen  with  cut-off  at  the 
position  of  9th  harmonic.  Noticeably  individual  harmonics  exhibit  tine  structure,  similarly  to  the  results 
reported  for  the  gas  phase  [106].  For  the  highest  intensity  30  TW/cm2  shown  in  Fig. 17(c),  a  second  plateau 
suddenly  emerges  with  cutoff  positioned  at  the  50th  harmonic.  The  appearance  of  second  plateau  with  har¬ 
monic  orders  extending  beyond  the  atomic  limit  of  the  corresponding  gas-phase  harmonics  indicates  strong 
dependence  of  HHG  on  the  band  structure  of  diamond.  Similar  results  are  found  in  HHG  spectra  from 
rare-gas  solids  [17]. 

Fig.  18(a-b)  shows  the  field-dependence  of  the  photon-energy  cutoff  for  the  first  and  second  plateau,  re¬ 
spectively.  In  contrast  to  the  cutoff  law  observed  for  gas-phase  harmonics,  cutoff  harmonics  generated  from 
diamond  bulk  display  almost  linear  scaling  with  the  field  strength  F.  This  linear  scaling  law  is  consistent 
with  experimental  observations  in  ZnO  [  6],  The  high-energy  cutoff  of  the  second  plateau  is  distinguish¬ 
able  for  strong  fields  with  F  ~  0.4V/A  corresponding  to  laser  intensity  I  >  10  TW/cm2.  The  slope  of  this 
linear  dependence  implies  that  the  cutoff  energy  increases  by  2  harmonic  orders  when  the  field  strength  is 
increased  by  0.5 V/nm.  In  contrast,  the  slope  of  the  photon-energy  cutoff  of  the  first  plateau  is  flatter  by  a 
factor  of  3.  Numerical  results  show  saturation  of  the  cutoff  of  the  first  plateau  for  field  strength  F  >  IV / A. 

In  Fig.  19  we  plot  the  field  dependence  of  the  spectral  intensity  of  individual  harmonics  Sn  from  the  first 
plateau.  The  laser  intensity  dependence  can  be  divided  into  3  regimes.  I.  Perturbative  regime:  The  3rd  and 
the  5th  harmonic  are  distinguishable  and  scaling  law  Sn  ~  Fn  is  exhibited.  II.  Intermediate  regime  (with 
F  ~  1  V/nm):  The  HHG  spectrum  extends  to  higher  orders  as  the  field  strength  increases.  The  intensity  of 
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different  harmonic  orders  show  the  same  power  law  Sn  ~  F(i  regardless  on  the  harmonic  order.  III.  High 
field  strength  with  F  >  5  V/nm:  the  intensity  of  individual  harmonics  become  indistinguishable. 

To  gain  insight  into  the  underlying  microscopic  mechanisms  of  hhg  in  diamond  we  plot  in  Fig.20(a-c) 
presents  the  time  evolution  of  the  current  density  distribution  for  three  different  laser  intensities  I  =2, 10  and 
30  TW/cnr.  In  the  perutbrative  regime  shown  in  Fig.20(a),  the  current  density  follows  quasi-adiabatically 
the  laser  vector  potential,  i.e. 


At) 


Jad(t)  =  ■  vk£„(k) 

v,k 


k=k(£) 


(39) 


where  k(f)  =  k  +  A(t)  is  shifted  crystal  momentum  and  e,.(k)  are  the  dispersion  curves  for  valence 
electrons.  For  the  intermediate  regime  in  Fig.20(b),  the  time  evolution  of  photoexcited  carriers  near  the  BZ 
edges  remains  adiabatic  due  to  the  wide  band  gaps  of  diamond,  however  states  near  the  BZ  center  do  not 
have  enough  time  to  fully  equilibrate  against  the  laser  vector  potential  A(t),  such  that  the  current  density  is 
strongly  distorted  after  the  peak  of  the  pulse.  This  distortion  expresses  the  non-adiabatic  effects  of  electron- 
hole  pair  excitation.  The  nonlinear  polarization  current  undergoes  rapid  changes  on  a  sub-femtosecond  time 
scale  tracing  out  the  burst  of  emitted  high  harmonics  after  the  peak  of  the  pulse.  For  high  field  strength  with 
intensity  I  =  30  TW/cnr,  the  prompt  excitation  of  electron-hole  pairs  becomes  highly  likely  as  shown  in 
Fig.20(c)  and  the  non-adiabatic  distortion  of  the  current  density  extends  over  higher  momentum  states  in 
the  BZ. 

The  Fourier  transform  of  the  even  part  of  the  current  density  /(k,  uj)  =  |  Jkw  +  J_kAJ  |  along  the  F-L 
line  is  shown  in  Fig.21(a-c)  .  The  odd  harmonic  structure  of  HHG  is  due  to  the  destructive  interference  of 
even  order  harmonics.  For  the  lowest  intensity  in  Fig.21(a),  sub-bandgap  harmonics  extend  from  the  BZ 
center  towards  the  edges;  the  recombination  of  electron-hole  pairs  near  the  BZ  center  gives  rise  to  intense 
emission  of  5th  harmonic,  while  the  7th  order  is  highly  suppressed.  For  higher  laser  intensity  of  10  TW/cm2 
in  Fig.21(b),  it  is  seen  that  the  HHG  spectrum  is  characterized  by  a  plateau  region  above  the  bandgap  with 
the  plateau  terminating  abruptly  at  the  position  of  the  9th  harmonic.  The  numerical  results  reveal  that  the 
plateau  and  the  the  high-energy  cutoff  converge  when  the  lowest  number  of  conduction  bands  included  in 
the  calculation  is  8.  Thus  we  confirm  that  HHG  depends  sensitively  on  details  of  the  band  structure  and  is 
characterized  by  strong  coupling  involving  multiple  single-particle  states.  For  the  highest  intensity  shown  in 
Fig.21(c),  the  spectrum  extends  to  high  harmonic  orders.  The  two  plateau  regions  can  be  distinguished,  the 
first  one  extends  to  15th  order,  decreases  gradually  in  spectral  intensity  with  the  increase  of  the  harmonic 
order  and  crosses  over  into  the  second  plateau  with  harmonics  extending  above  the  40th  order. 

To  non-adiabatic  electron  dynamics  associated  to  HHG  in  the  intermediate  regime  corresponding  to 
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laser  intensity  10  TW/cnr  is  shown  in  Fig.  22(a-c).  Fig.  22(a)  shows  the  temporal  profile  of  the  nonlinear 
current  corresponding  to  5-photon  transition  from  the  light-hole  into  the  lowest  conduction  band  (nc  =  1). 
The  current  waveform  follows  adiabatically  the  laser  vector  potential  during  the  first  paid  of  the  pulse  and  is 
distorted  2  fs  prior  to  the  peak  of  the  pulse;  the  subsequent  rapid  sub-cycle  oscillations  represent  the  build 
up  of  coherent  superposition  of  populations  in  the  light-hole  and  the  conduction  band;  steady-state  dipole 
oscillations  remain  after  the  end  of  the  pulse.  In  a  similar  way  the  laser-induced  hybridization  of  the  light- 
hole  with  a  higher-lying  conduction  band  (nc  =  2)  creates  rapid  dipole  oscillations  on  an  attosecond  time 
scale,  cf.  Fig.  22(b).  After  the  peak  of  the  pulse,  the  enhanced  steady-state  non-adiabatic  current  gives  rise 
to  emission  of  intense  7th  harmonic.  Fig.  22(c)  shows  the  nonlinear  current  corresponding  to  transitions 
into  and  out  of  the  nc  =  5  conduction  band,  noticeably  the  peak  of  the  current  is  delayed  by  5  fs  after  the 
peak  of  the  pulse,  the  time  profile  shows  that  the  9th  cutoff  harmonic  is  highly  virtual  and  exists  only  during 
the  second  half  of  the  pulse.  Similar  results  are  found  in  the  strong-field  regime  with  high  harmonic  orders 
due  to  transitions  into  and  out  of  highly  excited  conduction  band  states.  In  order  to  gain  further  insight  into 
the  microscopic  mechanism  of  the  harmonic  generation  we  separate  the  electron  dynamics  into  interand 
intraband  processes.  For  that  reason  we  expand  the  current  Eq.35  into  intraband  contribution: 


Jintra  (t)  =  ^2fnu{t)(e-vn  k)  (40) 

nk 

where  vn^  =  <9enk/<9(k)  and  interband  contribution  of  the  induced  transient  photocurrent 


J(M  =  X  ’  Prm'(k)  A(t)5nn ')  (41) 

nnf 

where  pnn/  is  the  matrix  representation  of  the  momentum  operator  and 

M 

Pnn'k{t)  =  y^{nk\'ipm(t))('ipm(t)\nfk)  (42) 

m= 1 

is  the  density  matrix  of  Bloch  electrons. 

In  Fig. 23  the  harmonic  spectrum  for  intraband  and  interband  currents  for  three  different  laser  intensities 
/  =2,10  and  30  TW/cm2  is  shown.  In  Fig.23  (a)  is  seen  that  for  the  subbandgap  harmonics  (first  and 
third  harmonics)  the  contributions  from  intraband  and  interband  currents  nearly  cancel  out  resulting  in  low 
total  current  causing  the  intensities  of  the  fundamental  and  the  third  harmonics  to  be  significantly  reduced. 
As  can  be  clearly  seen  from  Fig.23  (a)  (b)  and  (c)  the  higher  harmonics  (starting  from  the  fifth  harmonic) 
arc  due  to  interband  transitions  since  the  the  interband  current  predominates  over  the  intraband  current.  It 
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clearly  seen  that  the  spectral  intensity  of  the  intraband  current  decreases  drastically  after  the  21  harmonic 
order  (at  the  first  plateau).  At  the  second  plateau  its  spectral  intensity  drops  off  a  few  orders  of  magnitude 
in  respect  to  the  contribution  from  the  interband  current  as  seen  from  Fig.23  (c).  Therefore  we  conclude 
that  the  interband  process  is  the  dominant  one  for  the  harmonic  generation  after  the  5th  harmonic  order. 
The  time-profile  of  the  polarization  current  generating  the  5th  harmonic  for  the  three  different  intensities 
-  for  three  different  laser  intensities  I  =2,10  and  30  TW/cnr  is  shown  in  Fig.  24.  Clearly  harmonic 
generation  occurs  after  the  peak  of  the  driving  pulse,  the  pulse  duration  of  that  harmonic  is  shorter  than  the 
fundamental  period  T£.  The  fine  structure  of  individual  harmonics  is  imprinted  in  the  transient  modulation 
of  the  non-linear  electronic  current  on  a  femotsecond  timescale. 


D.  Conclusion 

We  presented  calculation  of  HF1G  spectrum  from  bulk  diamond  induced  by  an  intense  femtosecond  laser 
pulse.  Our  result  describes  the  intense  HHG  extending  from  the  perturbative  to  highly  non-perturbative  inci¬ 
dent  laser  intensity  regimes.  The  non-perutrbative  regime  of  HHG  occurs  for  laser  intensity  I  >  1 TW / cm2 . 
We  find  that  the  HHG  spectra  from  diamond  exhibit  two  plateaus,  the  second  plateau  extending  beyond  the 
limit  of  the  corresponding  gas-phase  harmonics  emitted  for  the  same  laser  intensity.  We  find  that  HHG 
and  cutoff  law  arc  sensitive  to  details  of  the  substrate  band  structure  involving  strong  interband  couplings 
between  multiple  single-particle  states.  The  numerical  results  make  evident  the  feasibility  of  diamond  for 
sub-femtosecond  pulse  generation. 
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L  r  x 


FIG.  1.  Static  band  structure  of  diamond  along  the  and  TL  directions  in  the  Brillouin  zone.  The  crystal  momen¬ 
tum  is  given  in  of  units  of  27r/ao- 


FIG.  2.  Time -evolution  of  the  electron  density  for  three  different  peak  intensities  of  30fs  laser  pulse  interacting  with 
diamond.  /  =  1.5  TW/cm2  (dotted  line),  I  =  10  TW/cm2  (dashed  line)  and  I  =  50  TW/cm2  (solid  line). 
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FIG.  3.  Time-evolution  of  the  macroscopic  polarization  density  along  the  laser  polarization  direction  for  three  differ¬ 
ent  peak  intensities  of  30fs  laser  pulse  interacting  with  diamond.  /  =  1.5  TW/cm2  (green  line),  I  =  10  TW/cm2  (red 
line)  and  I  =  50  TW/cm2  (black  line). 


FIG.  4.  Density  of  conduction  band  states  after  the  irradiation  of  diamond  with  30fs  laser  pulse  linearly  polarized 
parallel  to  bonds  between  carbon  atoms.  In  Fig. (a)  the  peak  intensity  of  the  light  pulse  is  I  =  1.5  TW/cm2.  In  Fig.(b) 
and  (c)  the  peak  intensity  is  I  =  10  TW/cm2  and  /  =  50  TW/cm2,  respectively.  The  kinetic  energy  of  the  released 
electrons  is  measured  relative  to  the  conduction  band  minimum. 
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FIG.  5.  Intensity  dependence  of  the  photoionization  yield  after  the  irradiation  of  diamond  with  30fs  laser  pulse 
linearly  polarized  parallel  to  bonds  between  carbon  atoms. 
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FIG.  6.  Perturbative  vs.  non-perturbative  multiphoton  ionization  of  diamond  interacting  with  30fs  laser  pulse  with 
intensity  =  1.5  TW/cm2.  Fig.  (a)  Time  evolution  of  the  modulus  of  the  Fourier  amplitudes  |ank|  in  the  expansion  of 
the  photoelectron  wave -packet  over  Bloch  states  with  k  =  (0,  0,  0)  by  neglecting  laser-induced  interactions  among 
valence  bands  (cf.  also  text).  Valence  and  conduction  band  amplitudes  are  given  by  blue  and  red  lines,  respectively. 
Fig.(b)  shows  the  time  evolution  of  the  modulus  of  the  Fourier  amplitudes  by  including  laser-induced  dipole-dipole 
interactions  between  the  three  valence  bands  at  the  Brillouin  zone  centre  (blue,green  and  black  lines),  the  photoioniza¬ 
tion  amplitude  is  given  by  the  red  line.  Fig.(c)  shows  the  imaginary  part  of  the  interband  polarization  by  neglecting  (red 
curve)  and  by  including  indirect  photoionization  pathways  (black  curve)  and  Fig.(d)  gives  the  frequency-dependence 
of  the  modulus  of  the  Fourier  transform  of  the  interband  polarization. 
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FIG.  7.  Comparison  of  theory  and  data.  Square  markers:  generated  charge  as  a  function  of  the  pulse  energy  from 
experimental  data.  Circles:  generated  charge  derived  from  the  theoretical  electron  yield  in  Fig.  5  distributed  in  an 
effective  ellipsoidal  volume  of  axis  wo  =  2.5  pm  and  zr  =  35  pm.  Black,  gray  and  light  gray  lines:  charge  calculated 
from  Eq.  15  as  a  function  of  the  pulse  energy  calculated  with  Eq.  15  for  wq=2,  2.5  and  3  pm. 


FIG.  8.  Time  evolution  of  the  free -electron  density  in  diamond  irradiated  by  200fs  laser  pulse  with  intensity  30 
TW/cm2,  linearly  polarized  along  the  (1,1,1)  direction.  The  red  curve  shows  the  cycle-averaged  electron  density  and 
the  blue  curve  is  the  electron  density,  (b)  The  cycle-averaged  carrier  densities  for  intensity  1=30,  40  and  50  TW/cm2 
are  shown  by  the  dashed-dotted,  dashed  and  solid  lines,  respectively.  The  position  of  the  pulse  peak  is  indicated  by  the 
vertical  dashed  line  and  Fig.  (c)  presents  the  number  density  of  non-adiabatically  excited  carriers  for  intensity  1=30, 
40  and  50  TW/cm2 
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FIG.  9.  (a)  Instantaneous  excitation  energy  of  electrons  interacting  with  200fs  laser  pulse  with  intensity  30  TW/cm2 
(linearly  polarized  along  the  (1,1,1)  direction),  the  green  curve  shows  the  cycle-averaged  electronic  excitation  energy, 
(b)  The  cycle-averaged  excitation  energy  for  laser  intensity  1=30,  40  and  50  TW/cm2  is  shown  by  the  dashed-dotted, 
dashed  and  solid  lines,  respectively.  The  position  of  the  pulse  peak  is  indicated  by  the  vertical  dashed  line. 


FIG.  10.  Density  of  conduction  states  after  the  irradiation  of  bulk  diamond  with  200fs  laser  pulse.  The  energy  is 
measured  relative  to  the  conduction  band  minimum.  The  laser  intensity  is  /  =10,  20  and  30  TW/cm2in  Fig.  (a-c), 
respectively.  The  vertical  dashed  line  in  Fig.a  indicates  the  threshold  for  impact  ionization 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


38 


„  1.2 

& 

CO 

c  0.6 

o 

CT1 

o 

“  0.0 

"c 

T3 

-0.6 

200  300  400  500  600  700  800 

time[fs] 


FIG.  11.  Fig.  (a)  Time-dependent  rates  including  carrier  generation  (positive  part)  and  laser-induced  recombination 
(negative  part).  The  laser  intensity  is  10  TW/cm2(dotted),  20  TW/cm2  (dashed)  and  30  TW/cm2  (solid  line).  Fig.(b) 
shows  the  intensity-dependent  impact  ionization  rate. 


time[fs] 


FIG.  12.  Conduction  electron  density  due  to  photoionization  only  (dashed  line)  and  including  the  impact  ionization 
(solid  line).  The  laser  intensity  is  18  TW/cm2.  The  vertical  dotted  line  indicates  the  position  of  the  pulse  peak. 
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FIG.  13.  Frequency  dependence  of  the  real  (solid  line)  and  imaginary  part  (dashed  line)  of  the  inverse  dielectric 
function  of  photoexcited  carriers  subjected  to  200  fs  laser  pulse  with  intensity  18  TW/cm2.  The  time  interval  is 
measured  relative  to  the  peak  of  the  pulse  ( t  =  0).  In  Fig.  (b)  t  -  30  fs,  and  in  Fig.  (c)  t  =  250  fs.  The  photon  energy 
is  indicated  by  the  vertical  dotted  line. 
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FIG.  14.  2D  intensity  map  of  the  logarithm  of  the  transient  efficacy  factor  of  laser-irradiated  diamond,  as  a  function 
of  the  normalized  (to  the  laser  wavelength)  LIPSS  wave  vector  components  (kx,  kv),  for  Instantaneous  bulk  plasma 
frequency  (a)  wp(t)  =  1.5 ujl  and  (b)  uip(t)  =  1.43wl.  The  laser  beam  is  linearly  polarized  along  the  y-axis  and  is 
normally  incident  to  the  surface. 
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FIG.  15.  Transient  variation  of  the  efficacy  factor  in  the  direction  perpendicular  to  the  laser  polarization.  In  Fig.  (a), 
the  dashed-dotted  line  the  EHP  frequency  u>p(t)  =  1.4wl  is  just  below  the  SPR  excitation  threshold,  the  solid  line 
represents  resonant  excitation  of  surface  plasmons  corresponding  to  instantaneous  bulk  plasma  frequency  ujp(t)  = 
1.43wl,  the  transient  increase  of  the  free-carrier  density  with  ujp(t)  =  1.45 u>l  (dashed  line)  and  uip(t)  =  1A7u>l 
(dotted  line)  results  in  red-shift  of  the  surface  plasmon  peak  towards  the  light  line  and  formation  of  near  wavelength 
ripples.  Fig.(b)  demonstrates  the  dependence  of  the  efficacy  factor  on  the  Drude  damping  time  re  =  100  fs  (solid 
line)  and  re  =  10  fs  (dashed  line).  The  laser  beam  is  linearly  polarized  along  the  y-axis  and  is  normally  incident  to 
the  surface,  A  and  A  designate  the  laser  wavelength  and  the  LIPSS  period,  respectively. 
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FIG.  16.  SEM  images  of  ablation  crater  edge  (ABL),  fine  (FR)  and  coarse  (CR)  rippled  regions  within  the  craters  on 
the  diamond  surface  for  N  =  10  (a),  30  (b),  100  (c),  and  300  (d)  pulses.  The  scale  bars  are  different  on  each  image 
and  the  bi-lateral  arrow  in  a)  shows  the  laser  polarization. 
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FIG.  17.  HHG  spectra  obtained  for  photon  energy  1.55  eV  and  pulse  duration  15/s.  The  intensity  is  2,10  and  30 
TW /cm2 
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FIG.  18.  High-energy  cutoff  as  a  function  of  the  driving  laser  field  for  the  first  plateau  (a)  and  the  second  plateau  in 

(b) 
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FIG.  19.  Dependence  of  the  spectral  intensity  of  individual  harmonics  on  the  electric  field  strength 
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FIG.  20.  Current  density  distribution  in  diamond  bulk  driven  by  intense  15fs  laser  pulse.  The  laser  intensity  is:  (a) 
I  =  2TW/cm 2,  (b)  I  =  10 TW/cm2  and  (c)  I  =  30 TW/cm2 
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FIG.  21.  Momentum-dependent  HHG  spectrum 
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FIG.  22.  Time  evolution  of  the  optically-induced  polarization  current  corresponding  to  recombination  of  electrons 
with  the  light-hole  band  close  to  the  BZ  center.  Fig.  (a)  5-photon  transition  into  the  lowest  conduction  band  (nc  =  1), 
Fig.(b)  7-photon  transition  into  the  first  excited  conduction  band  (nc  =  2).  Fig.(c)  presents  the  temporal  profile  of  the 
current  associated  with  the  emission  of  9th  cutoff  harmonic.  In  Fig.  (a-c)  the  laser  intensity  is  /  =  10  TW/cm2  and 
the  crystal  momentum  is  k  =  (0.007,  0.007,  0.007),  where  ccq  is  the  bulk  lattice  constant. 
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FIG.  23.  Harmonic  spectrum  for  interband  (solid  line),  intraband  (dashed  line)  and  total  (dotted  line)  currents 
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FIG.  24.  The  time-profile  of  the  polarization  current  generating  the  5th  harmonic  for  the  three  different  intensities  - 
for  three  different  laser  intensities  I  =2,10  and  30  TW/cm2 


Symbols,  Abbreviations,  and  Acronyms  are  explained  in  the  text 
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